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ABSTRACT 

This booklet is the ninth in a series of nine from 
the Teacher Training Institute at Hofstra University (New York) and 
describes the content and the approach of an institute course in 
which the participants use the personal computer as a personal tool 
within the mathematical discovery process of making conjectures, 
testing those conjectures, and verifying results and/or retesting. 
Pedagogical commentary and appropriate Pascal programs are presented 
for four topics, including: (1) iterated quadratic maps; (2) fractal 
geometry with applications; (3) conditional probability as applied to 
baseball league trends and results; and (4) recursion and induction 
in sorting routines. (JJK) 
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This booklet Is the last In a series of nine booklets which 
constitute the Hofstra University Teacher Training Institute 
(TTI) packet. The Institute was a National Science Foundation 
supported three-year prograa for exemplary secondary school 
■atheiatlcs teachers. Its purpose was to broaden and update the 
backgrounds Its participants with courses and special events and 
to train and support thee In preparing and delivering 
dissemination activities among their peers so that the 
Institute's effects would be multiplied. 

This packet of booklets describes the goals, development, 
structure, content, successes and failures of the Institute. We 
expect It to be of Interest and use to mathematics educators 
preparing their own teacher training programs and to teachers 
and students of mathematics exploring the many content areas 
described. 

"The Computer as an Experimental Tool In Teaching 
Mathematics" was a basic course offered as part of TTI's cycle 
of courses. This booklet describes the content and approach of 
this course - matfematlcs Is taught In a way analagous to how 
science Is generally taught. Using the personal computer as an 
experimental tool In a mathematics laboratory provides the 
student with an opportunity to get Involved In the discovery 
process: to make conjectures, to test them, to see the results 
and thus be able to adjust the conjectures being tested. A list 

o 

of course topics Is presented and several are described more 
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Science in high school and college is generally taught with the aid of 
laboratories. Laboratories provide the student with a chance to get involved in the 
discovery process: to make conjectures, to test them, and to see the results of this 
process. The advent of personal computers with graphics has made it possible to use 
the computer as a similar experimental tool in a mathematics laboratory. This paper 
describes the authors' collective experience in using the computer as an experimental 
tool in such a laboratory. 

This paper is an expanded version of a talk by one of the authors (HMH) at the 
December 1987 meeting of the National Council of Teachers of Mathematics. The 
paper reflects the experience of HMH in teaching an experimental course to selected 
high school teachers in a National Science Foundation Sponsored Teacher Training 
Institute at Hofstra University, as well as the experience of JB and RS in applying the 
techniques of that course in the high schools. We have not sought to provide a text on 
education, but rather to share some of our experiences. 

We thank Marie Hermann and Helene Morris for preparing this manuscript. 

Contents. 

1. Overview 

2. Iterated quadratic maps (HMH) 
a More on fractals (RS) 

4. Is there momentum in baseball (HMH and 1987 class) 

5. Sorting with the computer (JB) 
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1. Overview 

This set of notes is based on 1987 and 1988 class in the NSF-sponsored 
Teacher Training Institute at Hofstra University. The 1987 class met for a total of 30 
hours (twelve 2-1/2 hour classes); the 1 988 class for 9 such sessions. The audience 
consisted of superior, well-motivated high school teachers, who had a prior course in 
Pascal. Many had no recent experience with calculus. Many possible courses can be 
designed around these topics, which provide enough material and references for a 
one-year sequence. 

The goal of the class was to develop the use of the computer (in high school 
and calculus level mathematics) as an experimental tool in discovering mathematical 
ideas. The course emphasized experimental mathematics, in analogy with typical 
physics and chemistry classes. Thus we seek to use the "discovery method" in a 
variety of "advanced" topics accessable to high school students. Specific topics were 
chosen to demonstrate numerical and graphical techniques. The topics were chosen 
on the basis of mathematics level required, accessibility to the discovery method, and 
my personal interests, as well as to provide a useful and diverse experience for the 
audience. No specific attempt was made to cover the Advanced Placement syllabus or 
other syllabi, although the topic on sorting was added to the course at student request. 
A typical one-semester course would cover 4-6 topics. The topics are: 

1. A model for population dynamics. 

A simple model for population dynamics: the discrete time logistic model. Construction of the 
model. Method of computer simulation. Chaotic dynamics - determinism versus randomness. 

2. Random numbers. 

Linear congruence random number generator. What is a random number? Use of computer 
graphics. The birthday problem and its consequences. Probability and baseball - see Section 4. 
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3. Fractals I. 

Self-similarity. Regular fractals and recursion. See Section 3. 

4. Fractals II 

Fractals constructed by iterating quadratic maps. The Mandelbrot set and Julia set. Programming 
for time efficiency. See Section 2. 

5. Numerical integration. 

From a classical Greek formula for volume to Simpson's rule. Project on error analysis. 

6. Nun :al differentiation. 

Computer evaluation of limits, choosing the denominator appropriately. Project on error analysis. 

7. Hooke's law and the vibrating spring. 

Derivation of sinusoidal motion from elementary principles. 

8. Matrix models (age structured populations). 

Leslie matrices. Eigenvalues, eigenvectors, and convergence. Numerical computation of 
dominant eigenvalues. 

9. Sorting. 

See Section 5. 

We now illustrate several of the above topics. 
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2. Iterated quadratic maps 

The study of iterated quadratic maps is a nice demonstration of the use of 
computer graphics to uncover beautiful and interesting phenomena in mathematics. 

We begin by considering the process of iteration. Let f be a function from the 
real numbers to the real numbers. One may choose a real number xq , and form the 
sequence 

*0,xj «f <X0),X2=f (x 1 ),x 3 =f{x 2 ) 

(Many students have seen the iteration process in the computation of compound 
interest. Here the amount of money in an account at the end of each period is obtained 
by multiplying the amount at the end of the previous period by the quantity (1 + the 
interest rate per period). Thus compound interest involves iterating the function 

f (x) - (1 + r ) x 

where r denotes the interest rate and x the amount of money.) We are concerned 
with properties of the sequence of points obtained by iterating the map f. Among the 
simplest questions one can ask is whether the sequence xq,xi ,X2, ... is bounded or 
not. 

We also make this question simpler by first considering only linear maps: 
f(x) = ax+b. 

In this case, it is not hard to see that the sequence of points obtained by iteration is 

ERIC J 
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boundodif|a | < 1 , and in general unbounded if |a |>1. If b = 0, then f takes the 
simpler form f(x) * ax , and an easy calculation yields 

x n = a n x Q . 

The convergence properties of the sequence of points { x n ) follow easily. (In the 

example of compound interest above, a * 1 + r .) In the case where b is not 
necessarily equal to 0, 

x n = a n x Q + a n -lb + a n -2b + a n -3b + ... +b 

= a n x Q + (l-a")fc/ (1-a), 

summing a geometric series. The convergence properties of the sequence of points 
{ x n ) follow easily. We leave details and a more precise discussion to the reader. 

We now consider the case of quadratic maps: 
f(x) = ax 2 +bx +c. 

The answer is now more interesting, even in seemingly trivial cases such as f(x) = x2. 
Here the sequence of points obtained by iteration is bounded provided the first point 
x o satisfies | x 0 1 $ 1 , and unbounded provided that the first point x q satisfies 
I x o | > 1 . The most frequently studied case is 

f(x) = x 2 +c. 
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Here the results depend upon both the starting point x q and the parameter c 
in a complex way. 

We now consider iterating the equation 

f(z) = z2+c (2.1) 

for complex z (and possibly complex c). 

(We recall here that a complex number is a number of the form a + bi , where 
/ denotes V-1 . Complex numbers are added as if they were binomials with / as a 
variable; for example, (a + bi) + (c + di) * (a + b)+ (c + d)i . Similarly, complex 
numbers are multiplied as if they were binomials with / as a variable, except that / 2 
is replaced by -1 ; for example, (a + bi)(c+di) =(ac-bd) + (ad+bc)i. Complex 
numbers may be represented as points in the plane, with the complex number a + bi 
plotted as the point (a,b). The length of a complex number is then given by the 
Pythagorean theorem: \a + bi | » (a 2 + b 2 )™.) 

As above, for the results of iterating (2.1 ) for c - 0 are easy to see: if | z \ < 1 
then the iterates approach 0; If | z \ - 1 , then the iterates all have absolute value 1 , 
and if | z j > 1 then the iterates approach 00. The question is how do we study 
equation (2.1) if 0 is not zero. It is here that computers come to the rescue. 

First, to a computer a complex number is just a pair of real numbers or a vector, 

eg, 

z m x + iy corresponds to (x,y ). 
If we suppose that c is real, we may write the results of applying f to z ■ x + // as 
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Ynew = 2xy. 

Therefore one can just write a brief computer program to study the iteration, for 
example: 

program iterate; 

I^Xrjfihti^IJIarold M. Hastings. All rights reservei! except fliar 
C-Qoiesjnav be made for non-profit education! purples H\ a 

copy of this copyright notice appears on all such copies ) 

uses 

crt; 

var 

x,y,xnew,ynew,c : real; 
Umax : integer; 

begin 

writolnOHow many iterates ?'); 
readln(imax); 

writeIn(*What are the coordinates of the starting point, zO ?); 
readln(x,y); 
for i = 1 to imax do 
begin 

xnew := x*x - y*y +c; 
ynew := 2xy; 
x :-= xnew; 

1 2 
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y := ynew; 



writeln(i/ \x/ \y); 



end; 



end. 



Perhaps one can tell whether the sequence of points generated by this program 
approaches <» or not. However, since we cannot compute the entire sequence of 
iterates, it would be nice to have a criterion for testing whether the sequence of points 
generated by this program approaches Tha following lemma provides such a 
criterion. 

2.2. Lemma. Suppose |c | < 1 and \z | £ 2. Then |z 2 + c|£|z| + 1. 
Proof. | z 2 + c | £ | z 2 | - 1 c | (by the triangle inequality) 



as required.// 

Lemma 2.2 implies that once the length of any iterate z n exceeds 2, then 
subsequent iterates march off to |z n+ ^| s 2 + k. The following figure illustrates the 
inequalities in the Lemma 2.2. 



*l* 2 |-1 
2 2|z|-1 



(since |c|<1) 
(since \z\t2) 



= |z| + (|z|-1) 
2|Z| + 1 



(since |zj^2), 
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We now combine Lemma 2.2 with computer graphics to draw a picture of which points 
approach ~. The program is written in Turbo Pascal 5 for the IBM PC using Hercules 
compatible color graphics, and is easily modified for other languages and computers. 



program mandel; 

f Copyright (C\ Harold M. H astings. All rights reserved except that 
copies mav be made fo r non-profit educational purposes provided that a 
copy of this copyright notice appears on all such copies ) 



uses 



crt, graph; 



var 



ij: integer; 
k: integer; 
col: word; 
x,y; real; 
xn: real; 



z: 
c: 



real; 
real; 



(screen coordinates } 

(counter) 

(color) 

("math" coordinates) 

(nxew; the variable ynew is not needed in this 
program) 

(square of length of (x,y)) ; 
(constant in quadratic map) 



begin 



readln(c); 



(c must satisfy -1 < c < 1 ; the r der may add a 
check on c if required) 
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initgraph(0,0,' *); {Turbo 5 comand for hercules color graphics on 

IBM PC, modify as appropriate} 

fori:=0to319do 
for j := 0 to 1 99 do { loop over 320 x 200 graphics screen } 
begin 

x:=(i-160)/128.0; 

y := (100-j)/64.0; {x and y scale factors, adjust as 

appropriate, note the use of "100-j" in order to make 
the y-axis point up as usual, the action usually takes 
place within -1 < x < 1 and -1 < y < 1 } 

k := 0; {initialize counter} 

z := x*x + y*y; {initialize length squared} 

while ((k < 10) and (z < 4.0)) do 

{the "k < 10" criterion prevents infinite loops and 
might be adjusted by the reader. The "z < 4" criterion 
detects points for which we know that future iterates 
will approach infinity; see Lemma 2.2, above. Some 
points take more than 10 interations to escape the 
disk ,, z<4".} 
begin 

{first replace x + iy by (x+iy) squared + c) 
xn := x*x - y*y + c; {real part} 

y := 2*x*y; {imaginary part} 

x := xn; 
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z := x*x + y*y; {compute length } 

k := k+1 ; {iterate counter} 

end; 

if k := 10 then {point did not escape} 
col:=0 

else 

col := k - 3*(k div 3) + 1 ; {point did escape; 

color := k mod 3 + 1 
indicates time until 
escape. Other formulas 
might be used.} 

putpixel(i,j,col); 
end; { loop through points } 
end. {program} 



We now illustrate the results of several runs of this program, 
page 14: run with c =0.56. 

eags_15: run with c =0.56, but with "k < 1 0" replaced by " k <20", 
"if/c = 10" replaced by "if k = 20", and blown up by 200 %. 



17 



Hastings, Bernstein, and Silverstone 



October 19, 1989 



15 




ERIC 



BEST COFy AVAILABLE 

19 



Hastings, Bernstein, and Silverstone 



16 



We make several observations and suggest several exercises. 

(1) If c is real, as in the program, the resulting picture is symmetrical about 
both the x and y axes. We leave the proof as an exercise, but give a few hints. First, 
replace z by -z, and observes that z 2 + c = (-z ) 2 + c . Thus the fates of the points z 
and -z are the same. This yields symmetry about the origin. For symmetry about the 
y-axis, replace y by -y , and compute the length of 2 2 + c = (x +iy )2 + c and (x -iy )2 + c . 
Then show symmetry about the x-axis, using geometry. 

(2) What happens if c is not real? 

(3) The figures drawn by the program mandel above appear self-similar in 
that that each blob corresponds to two smaller blobs, of roughly similar shapes. The 
term "self-similar will be defined in Section 3, below. Here is a sketched explanation 
of why the figures are self similar. Suppose that a given blob consists of points which 
leave the circle |z | < 2 in k iterations. Then there are points which map to this blob in 
one iteration, and thus leave the circle |z| < 2 in /c+1 iterations. Since the map f {z) 
= z 2 + c is usually two to one, there will be two blobs which leave the circle |z | < 2 in 
/c+1 iterations. See the following figure. 



leave in k + 1 
Iterations 



leave in 

k + 1 iters. 




leave in k iterations 
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For self-similarity, we invoke the fact that f has a derivative V (z ) = 2z. Thus f 
rotates and stretches small vectors Az by a locally constant amount: 

f (z + Az ) = f(z) + 2z- Az + (small error). 

Multiplication by 2z multiplies lengths by 1 2z | and rotates vectors through an angle 
arg (z ) with cos (arg (z )) = x and sin (arg {z)) = y where z = x + // . We refer the 
reader to any introductory text on complex variables. This implies that blobs are 
stretched and rotated by locally constant amounts. 

We encourage the reader to experiment further. For further reading we suggest 
the following. 

References 

Mandelbrot, B.B. 1977, 1982. Fractals: form, chance and dimension, The 
fractal geometry of nature, Freeman, San Francisco. 

Pietgen, H.O., and P.H. Richter, 1986. The beauty of fractals: images of 
complex dynamical systems, Springer- Verlag, New York and Berlin. 



Remarks, We began with some simple questions. We first observed that the 
computer could be used to conduct experiments in order to try to answer those 
questions. We then saw that one could prove a useful rule (the lemma above) in order 
to use the computer to answer questions about the properties of certain sequences. 
This combines inductive reasoning from the experiments with deductive reasoning 
used in mathematical proofs. Many mathematicians are familiar with this combination 
of techniques, yet current curricula provide little useful experience for the student. 
Finally, we used computer graphics (combining algebra and geometry) to illustrate the 
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answers to a mathematical question. The study of the properties of the answers leads 
to a new field, the field of fractals, which is explored further in the next chapter. 



FRACTAL GEOMETRY 

AND ITS APPLICATIONS 
IN THE MATH CLASS 



Fractal Geometry is one of the fastest growing mathematical disciplines. 

The term M fractal" was coined by B, Mandelbrot of I.B.M, in the 

1970' s to describe the geometry of chaos. There are many ways to describe 

what a fractal is. Some of these descriptions are: 

a) A fractal is a figure that has a FRACTIONAL dimension. 

A line has- dimension 1 

A square has dimension 2 

A cube has dimension 3 

A fractal can have dimension 1.2345- 

b) A fractal is a process of becoming rather than being 

c) A fractal can be a sel f-similar object. 

This means that at any level of magnification, any p -t of a fractal 
can look exactly like the initial view. Self-similarity is related 
to the process of RECURSION 

d) An object of great complexity. 

e) A geometry that accurately describes the real world 
< or any other world ) 

In this session, we will examine the following ideas: 

a) Randomness and chaos... The CHAOS GAME 

b) Self-similarity , The KOCH SNOUIFLAKE and its cousxns 

c) Dimension, both topological and fractal 

d) Applications of these ideas in the Math curriculum 
e> Computer generation of these shapes 



THE CHAOS GAME 

Start with the vertices of a triangle, say X lf X a and X a . 

Let P be any. point in the plane. The CHAOS GAME is played as 

f ol lows: 
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a) Choose any one of the three vertices at random, say Xi . 

b) Mark the MIDPOINT of the segment P*~Xi . call this point P. 

c) Go back to step (a) 
Question: NHAT is the resulting object, and what are its properties 

NHY does the object appear as it does? 

OA 



Fractal < -J- > geometry 
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The figure produced by the CHAOS GAME is shown 
at the right. Although the object was formed 
by a random process, it appears to have a 
definite form and structure. 

Obser vat ions: 





a) The object has the property of being SELF-SIMILAR If you examine 
any of the sub-triangles, you will notice that they are IDENTICAL 
to the original, except for SCALE 

This property of sel f-si mi 1 ar it v defines the figure as a FRACTAL 

b> We can calculate the area of the object. 

A - A/4 - AC3/16) - AO/64) - - AC1/4X3/4)" - ••• 

where A is the area of the entire triangle 

c) This concept leads to further exploration: 

i) What happens i f we went 1/3 the way from the vertex? 

ii) What happens if we went 1/3 the way from P ? 

iii) What if a square, or any other polygon were ur »d instead of 
the triangle? 

The figure formed is called the SIERPINSKI TRIANGLE. Cousins of the 
triangle are the SIERPINSKI CARPET, pictured below 'left, and the 
MENGER SPONGE, pictured below right. 



By changing the initial conditions, but by using the same process, 

more realistic images, such as the fern, pictured below, can be generated. 







A simple BASIC program to generate the Sierpinski triangle is: 

Program is in AP PLESOFT. . Adjust graphics for your own computer ) 

10 X(l) : Yd) = : XC2) = : Y<2) = s XO) = : YO) « 

< Choose your own coordinates for the triangle's vertices ) 

20 XP.= s YP = : REM choose your own starting point (XP.YP) 

30 HGR : HC0L0R ■ 7 

40 R s INTO * RND(l) + 1 ) 

50 XP * ( X(R) + XP ) / 2 : YP = ( Y<R) YP ) / 2 

60 HPL0T XP,YP 

70 GOTO 40 or 

Fractal < > geometry 



THE KOCH SH0I4F LAKE 



The KOCH SNOWFLAKE is one of a class of fractals in which a straight 
line segment is replaced with a polygonal line, called a generator. 

The SNOWLINE is constructed in the following manner: 

a) Begin with a line segment 



Divide the segment into 3 equal 
parts, and replace the middle 
third with two line segments, 



d) 



repeat (b) on each of the 
resulting line segments. 



qo to (c) 




We observe that at each stage, the length of line segment increases by a 
factor of 4/3 , hence at stage n, the length of the "curve" is 

L„ = L*(4/3)" 

where L is the length of the initial line segment. 



Clearly, we can see that LIM L„ = 

n -> oo 

in a finite span. 



, hence is an unbounded length 



To do this program in BASIC is quite complex, because to do it efficiently 
requires RECURSION, or the ability of a subroutine (procedure , function) 
to call itself up. Below are two programs to draw the Koch line and 
the snow flake: 



LOGO 

TO KOCH : LENGTH : LEVEL 

IF : LEVEL = 0 [FORWARD : LENGTH STOP] 
KOCH s LENGTH/3 : LEVEL- 1 
LEFT 60 

KOCH : LENGTH/3 : LEVEL- 1 
RIGHT 120 

KOCH i LENGTH/3 :LEVEL-1 
LEFT 60 

KOCH : LENGTH/3 : LEVEL- 1 

END 



PASCAL 

PROCEDURE KOCH <D: REAL; L: INTEGER) 
BEGIN 

IF L = 0 THEN MOVE(D) 
ELSE 
BEGIN 

K0CH(D/3,L-1);TURN(G0); 
KOCHCD/3, L-l ) ; TURN (-120) 
KOCHCD/3, L-l ) j TURN (60) ? 
K0CH(D/3,L-1) 
END 
END; 
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The Koch Line hats some interesting properties: 



a) For any initial span, the "length" of the curve is infinite 
b> Although the curve is continuous, at no point is there a derivative 
< an example of where continuity is not sufficient for the existenc 
of a derivative. ) 



The koch line can be expanded to the KOCH SNOWFLAKE by: 

LOGO PASCAL 

TO SNOWFLAKE : LENGTH : LEVEL FOR I := 1 TO 3 DO 

REPEAT 3 CKOCH : LEVEL : LENGTH RIGHT 1203 BEGIN 
END KOCH < LENGTH , LEVEL ) ; 

TURNC-120) 
END; 

The first five stages of the SNOWFLAKE are shown below: 



\ 



\ 



v 







Can we prove that the area is bounded, without actual computation? 
Can the "actual" area of the snowflake be calculated? 

Here we have an example of an infinite perimeter bounding a finite area. 

Fractals are used to model the real world. One application of the snowflak 
is in the representation of islands and clouds. The coastlines of islands 
and the boundries of clouds can be thought of as random. By introducing 
rnadorn lengths and angles into the generator, we can simulate these natura 
structures. 




Random KOCH ISLAND 
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DIMENSION 

It is our normal understanding that the "dimension" of an object is an 
integer. We know this because we asssociate dimension with direction. 
We say that a line is 1 -dim because we can travel in essentially one 
direction, forwards Cor backwards). A plane is 2-dim because we can walk 
not only forwards < & backwards) but up and down. We have a sense that 
dimension is related to the "amount of space that is taken up". This 
notion of dimension is called the TOPOLOGICAL definition of dimension. 

We can see how dimension is developed by the following considerations: 

a) consider the line of length 1 ______ 

Suppose that we triple the length. Now we can see that 
three of the original segments will cover the new length. 

We can write 3 l ■ 3, where the power, 1 , is the dimension. 

b) I f we consider a unit square, and triple the dimensions, 

it requires 9 of the original squares to cover the new one. 

We can write 3 a = 9, where the power, 2 , is the dimension. 

c) If we consider a unit cube, and triple the dimensions, 

it requires 27 of the original cubes to "fill" the new one. 

We can write 3 3 ■ 27, where the power, 3 , is the dimension. 

These examples lead to a relationship between the number of self-similar 
parts N generated and the scaling factor <R). It is: 

pat °" a N or R p = N or D - log (N) /loaCR) 







R 


N 


D 




Unit 


1 i ne 


3 


3 


1 


3 l = 3 


Unit 


square 


3 


9 


2 


3* ■ 9 


Unit 


cube 


3 


27 


_ 


3» - 27 



and these conform to our usual sense of dimension. 

Applying this relationship to our two "strange" objects, the gasket and tht 
snowflake, we see: 

TRIANGLE If we think of the lower left as the "unit" triangle, 
the triangle's sides are twice as large, and there are 
three triangles generated. Hence 

R = 2 N = 3 Log<3)/log<2> ■ 1.5849- 

SNOWFLAKE Each length is divided into 3 equal sections, and 4 segments 
replace the original 3, hence 

R = 3 N ■ 4 Log (4) /LogC3) ■ 1.2618- *• 

These strange objects have fractional dimensions. A question that arises 
is: "What is a meaning of a fractional dimension?" 

ERIC 

vMomam Fractal < -5- > geometry 



The CARPET'S dimension can be calculated by by observing that 

the original square is reduced by a factor of 3, and 8 squares are 

generated: 

r s 3 N = 8 log<8>/log<3) = 1.89- - 



THE CANTOR SET 

One of the earliest fractals developed was the CANTOR SET. It is derived by 

a) Start with a line segment 

b) Divide into 3 equal parts, and 

remove the middle third 

c> on each of the resulting segments, 
repeat part <b> 



FRACTAL GENERATION USING TRANSFORMATIONS OF THE PLANE 

The regular fractals discussed so far can be generated by considering 
the movements of points in the plane < or space ) by means of 

AFFINE TRANSFORMATIONS 

An AFFINE TRANSFORMATION is a LINEAR TRANSFORMATION followed by some 
TRANSLATION or. SHIFT . LINEAR FUNCTIONS consist of such movements as 

REFLECTION EXPANSION DILATION ROTATION SHEARING 

Suppose that T is a linea^ function of the plane to itself, and suppose 
that the point <x,y) is mapped, under T, to the point <x',y' >. 
We can then write: k 

<x,y)T *> <*' ,y' ) * 
<x',y') is called the IMAGE of <x,y> under T 

This means that trere are real numbers a, b, c and d, where 

x' * ax + by 
y' - cx + dy 

This system of linear equations can be written as a MATRIX equation 



<".x' ,y' ) = <x,y)T = <x,y)* 
An AFFINE TRANSFORMATIO N. A, can ue written 



(: :) 



<x,y)A = Cx,y)T + <:r,s> 

where Cr,s> is a trarslation, or shift r units Horizontally and s units 
vert ical ly . 

o 
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Returning to the SERIPINSKI TRIANGLE: 

Our goal is to see where points go under affine transformations. 

For the TRIANGLE, we will need to consider THREE transformations 

Ai : The image of (x,y) will be in the LOWER-LEFT triangle. 
Aa s The image of (x,y) will be in tf- > LOWER-RIGHT triangle 
Aa : The image of (x,y) will be in the TOP triangle. 




This qives the following set of affine transformations: 



For A 4 



(0,0) A* = (0,0) 

a,o)Ai = <i/2,o:> 

<h,k)A» = (h/2,k/2) 



IM0 




For Aa 



Resulting in: (x,y)Ai = (x,y)* 



(0,0)A a = (1/2,0) 
(1,0) A a = (1,0) 
(h,k)Aa = ((h+l)/2,k/2) 




A. 




For Aa 
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Resulting in: (x,y)A a 

(0,0)Aa = (h/2,k/2) 
(l,0)Aa ■ ((h+l)/2,k/2) 
(h,k)Aa - (h,k) 



fo,o) 



Resulting in: (x,y)Aa = (x,y)* 
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+ (1/2,0) 




0,o) 



A/2 o\ 

/ + (h/2, 

\0 1/2/ 



k/2) 



SOME APPLICATIONS OF FRACTALS 



Fractal analysis is being used in studying almost all natural phenomenon. 
One example is in Motion pictures. The STAR TREK AND STAR HARS worlds 
were generated by fractal programming of computers. Below is a "simple" 
example of a randomly-generated landscape. 




A GALLERY OF FRACTALS 



a) From the chaos game, if the point is chosen a distance 1/3 to the 

vertex instead of the l/2 f the result is: 

A 
A A 

A 

A A A A 

* A 

.% A A A 

A A A A 

A A AA A A A •« 



b) The world-famous BIFRUCATION graph 




•:> The grand-daddy of all fractals, th« MANDELBROT SET 




BEST COPY AVAILABLE 



A SET OF PROBLEMS FOR YOU TO TRY 



Start with a square and replace each side 



of the square with the generator 




Let P = initial perimeter and A - initial area. 

a) Letting SGKO) oe the initial square, draw SQC1) and SQ<2) 

b> Find the dimension of this fractal < SQ(— ) ) 

c) Find the area ACn) of SQ<n) and hence AF = Lim A<n) 

n->"> 

d) Find the perimeter PCn) of SQ<n) , and PF = Lim P(n) 

n-> «© 

e) Write a program to produce the this fractal. 
This fractal is known as the KOCH ISLAND 

Give a geometric argument to show that the area of the SNOWFLAKE 
is bounded by the hexagon that circumscribes it at level 1. 



Refer to the CANTOR SET, page 6, 

a) Calculate its dimension 

b) calculate its length 

Start with a square and form the fractal by replacing each of the 




sides 



with the generator 



a) Draw level 1 and level 2 of this fractal. 



b) What is the dimension of this fractal? 



Do the same as problem 3 with the generator 
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Refer to the SIERPINSKI CARPET on page 2 

Let MCO) be the complete square at level 0, 

M<1) be the square, at level 1, with the central square 
removed, 

M<2> be at level 2, where the 8 additional squares are 
removed. 

a) What is the dimension of this fractal? 

b) At. stage n 

i) SCn> = number of squares removed. What is S(n>? 
ii) What is the area of each of these squares ? 
iii) What is the total area, TA<n), removed 

c) What is the area of this fractal ? 

Refer the the MENGER Sponge on page 2. 

a) Find the fractal dimension of the sponge 

b) Calculate the surface area of this fractal. 

c) Calculate the volume of this fractal. 

Consider the following construction ( drawn in 3-D perspective) 
level 0 Level 1 level 2 




a) 


What 


i s 


the 


dimension of this fractal 




b) 


What 


is 


the 


surface area of this fractal? 




c :> 


What 


i s 


the 


volume of this fractal? 




d) 


What 


i s 


the 


height of this fractal? 




e) 


What 


i s 


the 


relation between this fractal 


and the MENGER 




SPONGE? 








a) 


Derive 


the 


affine transformations for the 


CARPET 


b) 


Der i ve 


the 


affine transformations for the 


SPONGE 
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SOLUTION SUGGESTIONS 



1. 



a) SGK1) is 




b) 

.; ) 

d) 

*:> 



D = Log(8)/log(4) * 1.5 
A(n) = A for every A. 
P(n) ~ 2t*P — > «o 
see next page (after # 9 



3. 
4. 



8. 



Just prove that at each stage the triangles added in are within th 
line joining the furthest endpoints of the two adjacent sides: 




a) D = Log (2 :> /Log (3) ^ . 6C 

Level 1: 

Dim = log (8) /log (3) 1.89 

Level 1: 



Dim = log (9) /Log (3) = 2 
( Is this a fractal???) 



b) 0 



a) D = log(8)/log(3) 

b) i) S(n) = 8" 

c) 0 



ii) (1/9")*A 



iii) TA(n) ■ (8/9)"*A 



( consider the geometric series described by part (b) 

a) Log(20)/log(3) Z 2.7268 

b) At each stage, each sub-cube's surface area is increased by 
a factor of 2 

c) 0... relate to CARPET 

a) log(13)/log(3) * 2.334 

b) At stage n, the surface-area is (13/9)" of the initial square 
hence is infinite. 

c) At stage n, the volume is increased by < 13/27) "V 27 
(assume there is is an initial cube of volume 1) 

Volume = 1/14 of the initial surrounding cube 

d) Height is (1/3) + (1/9) + (1/27) + ••• = 1/2 



ERIC 
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a) Requires 8 affine transformations 

IT 






s 


s 
















/ 



b) Requires 20 affine transformations. 



Note the geometric "similarities" between the following fractals. 
Do you suppose that there may be some relationhip between them? 






LOGO program 

TO ISLAND : LENGTH : LEVEL 

IF : LEVEL = 0 FD : LENGTH STOP 
ISLAND : LENGTH/4 :LEVEL-1 
LEFT 90 

ISLAND : LENGTH/4 : LEVEL- 1 
RIGHT 90 

ISLAND : LENGTH/4 : LEVEL- 1 
ISLAND : LENGTH/4 : LEVEL- 1 
LEFT 90 

ISLAND : LENGTH/4 : LEVEL- 1 
LEFT 90 

ISLAND : LENGTH/4 : LEVEL- 1 
RIGHT 90 

ISLAND s LENGTH/4 :LEVEL-1 

END 




KOCH ISLAND 60 2 



TO KOCH ISLAND : LENGTH : LEVEL 

REPEAT 4 C ISLAND : LENGTH : LEVEL RIGHT 90 ] 

END 
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ADDENDUM 
I 



n.e following problems are examples of how the idea of U£LF -SIMILARITY 
can be used to simplify the? solution. 

!• C ONTINUED FRAC ;;)NS 

Let X = 1 + I 

1 + I 

1 + 1 

1 + 



The sel f -simi larity in this problem is quite evident. We can 
re-write this expression as: 

X = I + I 
X 

which results in the quadratic equation y; 2 - X ~ 1 = 0 
hence X = 1 + /5 



This is a problem that was given to me by a student: 

Ln(X 

Ln<:x 

Evaluate Ln(X 



Let thi.-. expression be represented Ty Y, then clearly, we can write 

Y = Ln(X v > 

Y = Y-LruX) 
1 = Ln(X) 

e = X 
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• »■> .«1 I.jm; n .jo ; n"b i tii- w* 'C g ; "o hiv !jV Hi '"ic«''i !Jr -wn ■ f 

► ti; A- Kit oy . !" <.otf r f NH. 

A arj B, m turn, roll a fair Jie. A rolls first. i.'.e first persoi 
to oil a "6" wiria. What i«i the probability, a, that A wins? 

The? tree below shows the proqresD of the game for the first 5 rollj, 




Clearly, a * (1/6) + (25/36X1/6) + (25/36) 2 ( 1/6 > + ■ • ■ 

We can use the idea of SELF -SIMILARITy by observing that the 
tree can be re -writ ten to look like 



giving the equation a = (1/6) + <25/3b)a a =■ 6/11 




The classic GAMBLER 7 3 RUIN PROBL EM 

A and B play a game. A starts with $3 and B starts with *2. The 

game is as fol lows: 

Each, in turn, flips a fair coin. If the tosser gets a HEAD, 
then he receives '£1 from the other person; otherwise he gives 
£1 to the other person. 

What is the probability that the first person wins? 



Tn the solution, a tree will again be used. Let (a,b) represent 
how much money each has as each stage of the game. A tree, for the 
first 5 tosses is shown below: 




which gives the infinite series 

a = (1/4) + 2<1/16) + 5(1/64) + 13(1/256) + 



which is not an easy series to evaluate. 
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which gives the infinite series: 



a = (1/4) + U/4>a + (l/4)a + Cl/16)a + Cl/64)a + 

= (1/4) + (1/4) a + ( 1/3) a (Sum the infinite geometric series) 
Hence a = 3/5 



We can reduce the tree still further by making the observation that 
the state (2,3) is the same as the state (3,2), except that it 
represents the probability of B winning < b = 1 - a ) 




This now yives the equation a = (1/4) + (l/4!>a + (l/2)<l-a) 



The class of problems where we can use the techniques of SELF -S IM I LARITY 
is large and varied. This method of solution introduces the student to 
another way of viewing phenomena as a process, and in a fractal way. 
This method can also help to introduce the student to the ideas and 
techniques of RECURSION and and to reinforce the tool of INDUCTION* 
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I J 

: liov,aw i'i colltv>:on of AFFINE TRANSFORMATIONS .ha: you can 
in construction yc r v. -j 1 f -sifluliar fractals. 



Let '-.,.) be the given point, ar-d (x',y') = (x,y)A 



SCALING 



CD 



ROTATION 



A = 



/ cos(a) sin(a)\ 
\ -sin (a) cos (a) J 



REFLECTION 



a) X-axis 



b> Y-axis 



A = 



A = 



C .:) 

r. i) 



SHEAR 



a.' X -direct ion 



b> Y- direct ion 



Izz 



b 



TRANSLATION (note: write Cx f y as Cx,y,l) ) 



o 



A = 



A = 



A = 



f: ;) 

OH) 
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Is there momentum in baseball? 



Harold H . Hastings 
and 

The 1987 TTI Class* 

Department of Mathematics 
Hofstra University 
Hempstead, NY 11550 

Sportswriters often discuss "momentum". We sought to 
determine whether there is momentum in baseball , in one 
special case, the 1986 Mets , using elementary probability. 
The results are interesting and the problem and techniques 
interested our class (outstanding high school teachers 
participating in a National Science Foundation sponsored 
Teacher Training Institute ) and should interest many others . 
The calculations may stimulate additional investigation into 
ordinary events using elementary probability . 



The 1986 Hets season** may be summarized as follows, 
where Wn denotes a string of n Wins , and Ln a string of n 
losses 5 

W 2, L 3, W 11, L 1, W 6, L 1, W 1, L 2, W 1, L 2, W 3, L 2 
W 6, L 2, W 1, L i,W2,L 1, W 3, L 1, W7,L2, W 1 , L 1 
W 1, L 2, W 8, L 1 , W 1 , L 3, W 5, L 3, W 3, L 2, W 3, L 2 
W 3, L 2, W 3, L 1, W 3, L 2, W 1, L4, W 6, L 1, W 5, L 1 
W i.Li.W 4,La,W 1 , L 4, W 3, L 1,W i,L i,W2,L 1 
W 4, L 1, W 5 (108 wins in 162 games). 



We shall say define momentum as the tendency for wins to 
follow wins . This allows the questions about momentum to be 
phrased mathematically as follows : 



1 ) If the Mets won the last game , what is the probability 
that they win the present game? 

3) If the Mets won the last n games (for a fixed n) , what is 
the probability that they win the present game? 

3) Are there any statistically significant differences? 

We shall use test the null hypothesis that there is no 
momentum , in which case the results of successive games are 
independent . More precisely , we consider the Bernoulli trials 
(see [1]) or any elementary statistics textbook) model that 
the 1987 Mets season consisted of 16a independent games, with 
a probability p of winning each game of 108/163 or .667. 

F'or n independent Bernouilli trials , each with a 
probability p of success and q = i-p of failure, and np and 
nq i a5, the expected results are essentially normal with 
mean np and variance npq [ 1 ] . 

There were n = 107 games played following wins (the Mets 
won the last game). Of these, the Mets won 76, this gives a 
conditional probability for a win following a win of 

P i win j previous win) = 76/107 = .710. 

However , the null hypothesis gives an expected number of wins 
of 



np = 107 x .667 s 71.3, 



with a standard deviation of 

fnpq = /167 x .667"x .333 = 4.88 

games . The number of excess wins , 4.7, is less than one 
standard deviation , and thus not statistically significant . 

These calculations were repeated for the case of 2 or 3 
previous consecutive wins , with the following results . 



Number 
of pre- 
vious 
wins 

2 
3 



Wins 
(trials) 

75 
53 



Actual results 



(successes ) 

54 
35 



Prob. of 
win 

.720 
.660 



Null hypothesis 



Expected 

wins 

np 

50 

35. 3 



Standard 
deviation 
fnpq 

4.08 
3.43 



In neither case is the difference between: the null hypothesis 
and the actual results statistically significant . 

We invite the reader to continue with larger n , as we 
did, and see what happens. (We found no statistically 
signif ican "momentum" with our definition . ) We note that 
had we found momentum from one game to the next, we would 
have considered a MarKov chain model (see, for example [1]) 
which estimates the probability of winning each game from the 
results of the immediately previous game. We also invite 
other definitions of momentum . 
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3. SORTING WITH THE COMPUTER* Joyce Bernstein 

t r i * v odik'Z t i on t o Gor t :: 

^ 3 aht ?C) Jo t> Ser nst *i n f J 038 
T^achi,;*:; :v..rt .-.*lgo. *bhm;i i^. a :; ;.inpui >v aci&nc^ clc**^ ovicides a. 
nice briciqc* b»;t wo*n &v ;jq. ammi ng mathematics. It becomes necessary 

to understand what graphs of certain functions look like far removed 
from; the oriqin. Recursion, used in the most effecient sort, is & 
close cousin of induction 

When choosing a sort, there are three main consider at ionsn 
is Programming time - for small arrays or files, use a simple »ort« 
ii« Execution time - measure of effeciency, a function of number of 
comparisons and data movements* 

i i i s Memor y r equ i r ernent s - usuaJ 1 y sac v i f i c ed i n a t r ade~o f f for 
!j ^- L »: x>r" f f i' i enc y « 

~ff e-r.cy ht*<: - .v,* ::: Lmpov fcar.t a : : the nusTib^ / of -file-;* it-, bo *ort«ri 
ijjw-fc?.; 1c .*v.->m Ui: m-sr^^ure wrffeciency using <& fLy,u'iO p v-opo. i o; i^l to some 
fUii".t.;.on gCN> f of N f thr. -*i e of the sort. "B.i.ci 0"' not&ti^n, 
C Cg v i -1 -> ..' , mparts proportional to q (N) , 

fCNi = C(g<N>) «> |<>0 IfCNJ! < k! yCHJ! 
For all of the sorts we will consider, which require both comparisons 
and data movements, g(N) is either n* or niog a n„ Ti .^tractive 
to graph n a and nlog a n, especially as n gets large. 



I. Exchange sorts: Exchange sorts move data int.- place-, by position, 
one at a time. 



.y, .Sort f one of l.he easiest sorts to program, brings the 
smallest item to the front of the list in a manner which resembles 
"bubbl ing". 

Th e a 1 g or i t h m d o e s t h e f o 1 J o w i n g s 

is Compare items in pairs from the top of the list to the bottom, 
exchanging them if th* front item in larger than the back il«in. 
(. n - 1 c omp ar i son s ) » 

iia Repeat the process another n~2 times. The jth iteration makes 
n~J comparisons. We can sum the comparisons, n <n~i )/2, and assume 
that swaps occur about half the time, or n<n~-i)/4. 

!n(n-l)/4l * C 1/4) I n*-n ! £ (1/4)^*1 » 0<n«) 
Bubblesort is obviously an DCn*) sort."" 
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The Q yicfcsort . invented by C.A.R. Hoar e in 1962, moves one record at a 
time into its final position. However, in the process, which is 
recursive, the other records are moved to a position closer to their 
final onou. The algorithm involves a pardoning of the array into 2* 
:^ i o.v.r; or, th\? .jfch .:al I in- a i ■ In bhe first pas:>, a pi /...,t element is 

. - -., ...J (ui.:. a! : • t :.!■:<■ .m! of U.e fiiui this 

wlcM.i^i.l : ^ .i i. t .> ? L ."a! |jo.:*l bior'iy -all %-l ttftwotsi iw.».-jllev than it are 
:Jiiv-sc! of anr' a] 1 iiciriK. T .K-r *;*r^ l.-ihiru - ' it, 7h< <*] gor i thin c^I h 

£ f itself, uhoosing a pi. vol froui each "half" of the list, arid the simpler 

problem of sorting two smaller arrays is set up. The comparison takes 
place as follows! (assume the first pivot is item 1. When would you 
> want to make a different choice?). Assign two variables, i&j* with i 

initially 2 and j initially N. Compare element i to the pivot. If it 
is smaller, increase i by 1. Stop when you get to an element greater 
than the pivot- Do the reverse with item j. Keep decreasing j by 1 
until you get to an item less than the pivot. Exchange items i and j. 

The quicksort i:; a powerful □(nlog^n) a»or I . Drily necessary rr.oves are 
it'&d .' f v ' du '. \ i . s": o'r" ..- 1' !* ■ , fcw: 1 H ! 1 . t.'. c. 1 iijcv !i v!..v. can l*e I'liadr.'' ^v^n mcr* 
powerful by wvi.i^ pointers :l ynJ of data. 



II „ Insertion sorts: Each successive item is moved into an already 
sorted list. 



Linear insertion mimmioks the way many of us arrange playing cards in 
our hand when we pick fchem up one at a time. Pick up any card. Pick 
up the second card and put it in order- Pick up the third c<*.rd and 
put it in order, etc. The algorithm inserts the new element j by 
first comparing it bo element j - 1, which is already sorted, then to 
element j - 2, etc, moving each of these elements up one place until 
it f i n ds an e 3 e men t sma 1 1 e r !; h an i t . T* .h. men fc i s p 1 *u: ed i n its 

proper place. Notice the special condition which takes place when the 
element to br inserted i«; the smallest element on the Hot, i;p to the 
time oT thia insertion. The while loop is bypassed, so that the 
rompar L^on a!" place!] > tuinp i, not read w"hen place - 0. 

In order bo insert the .jth item, approximately CJ-D/2 coiup^r i ^on» and 
data movements are required. Summing over the n elements in the 
array, we see that linear insertion is an QCn'*) sort. It is a very 
effecient sort when the data is already almost in order. 

glLW.:/ ln£g£tj&Q is a simple, effecient sort. Each new element is 

inserted into the previously sorted partial list by using a binary 
search to locate the place of insertion. As with all "divide and 
conquer" strategies, binary insertion is an 0Cnlog a n) sort. 




91 m l Effo rt i<s a powerful insertion sort. Studies of N randomly 
arranged d^fca it^w:> ..ln-.w U.*-t v-2cm«nts y on the average, travel a 
distance of N/3 places to fi.iai, sorted positions. In 1959, Donald 
Shell used this fact in his: r or t algorithm, which initially moves data 
over larcje distanced, aptf- ^ ;imat ely equal to N/3. This method has the 
tendency, on the aver «*qe f nf moving elements closer to their final 
location very early in the -sort process. 

Shell partitioned the list to he sorted into k "chains" , each 
approximately 1/3 the <->ise of the total list. Each' chain was "sorted" 
using an exchange Chubb! The algorithm shown here, somewhat 

improved, uses a linear insertion for each chain. The first chain 
consists of elements 1, k-K, 2k+l, etc. The second chain consists of 
elements 2, k+2, 2k+2, *tc.» and the Jth chain <j<**k) consists of 
elements j, ,i*k, j*2k, etc. Thus, in the initial set of passes, 
elements are separated by a distance of k. K is then decremented to 
approximately k/3 *nd the proi/ess is repeated. When k is one, the 
list ^should be nearly irjir t * h -*.id a final linear insertion sort of the 
entire lint or.cuts. <*• linear y.ort at the end is faster than a binary 
sort bd..cvi<.;e v although binary m^tho.Js are better for randomly sorted 
larqe a linear nor fc ic- faster if the number of comparisons is 

smal I or i f the 1 L st is al most i n order (see note) . 

There are various versions of this sort, each using different rules 
for partitioning the list. A less effecient algorithm uses N/2, N/4, 
etc. The ant-* shown *hore is better than this version. Another 
version of the Shell sort uses decrements of -..,2^-1, • 
...,31,15,7,3,1. The number of compar isions, and moves is relatively 
low. Note that each p-*jr of cuccessivo decrements is relatively 
prime. 

note: It is interesting tu note that for small n, an CKn 3 *) sort is 
often more effecienct than an OCnlog^n) sort because their simplicity 
requires little ovc ' A^epl comparisons and passes, making the 

constant; of proportionality relatively small. 





BEST COPY AVAILABLE 




Program Sorts; ( input 9 output) J 
{Joyce Bf?rn3tc»i 1 il 
{F*b* £?P» 1?88> 

{Program which riemostrate«a an improved hiihblesartu quicksort,, 1 jne*r} 
{insertion sort., binary insertion sort* and shellsortl 
{copy for educational use only} 
const 

mexnum ** SOOO? 
type 

list « array Cl.«maxnum3 of integer? 
var 

n ? {size of array} 

choi ce« integer ? f menu selection! 

a jh » list? Cunsorted and sorted lisbs! 

{ * **##*##***#******procedure makel xst*#*###*#**#***^**********#******.J 
procedure makelisbtvar aslists nsinteqer); 

{generates an array of size n of random numbers greater than o< squall 

(zero and less than 1000} 

var 

i a integer s 
begin 

for i «« 1 to n do 

aCi3 j» random ( 1000) ? 

end S 

C**##4HMHt***#*##**#procedHre showlist**************^ 
procedure showl ist(a: 1 1st sin: integer ) ; 
{displays array in rows of 10 elements} 
var 

i 8 integer ; 
begin 

for i 1 « 1 to n do 

begin 

writ@<ati3n^ji * 51 > s 

if i mod 10 « 0 then writeln; 

end * 

wr i tel n? 
end 5 

{##***##**#*#*#*#*#prac®duv® insert ion*****************^ 

procedure i nser t ion (var a: 1 ist Sns integer ) 5 

var 

i ? {list index for move once place is found! 

.in {index for array item begin placed! 

place, {eventually, index of first item smaller than item to insert} 
temps integer? {value of item being inserted! 

founds boolean; {flag for location of insertion spot! 

begin 

for i •» S to n do 
begin 

temp 8« at J {first j ~ 1 elements already sorted! 

place 1* \ -J J 
found 1* false? 

while (place > 0 ) and not found do 
if atplace3 > temp 

then place i« place - 1 
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e I se found 8 « true? J 
for i i - t dewnto place + t do 

a r \ + n s * a r i :i j 

afplace 4- 1T ?^ temp 
end 
end J 

{*^^*#^#^*#^*#**prDceclur© binaryinsertion**************^ 

procedure binaryinsertion< var a; list? nsinfceqer); 

var 

{list index for move once place is* found > 
j? {index of item beinq inserted} 

topjbottomjimiddle^ {section boundaries and center} 
temp « i nteqer $ 
beqi n 

for j »=* J? to n do 
taeq i n 

temp ?*= aT j 1 ; 
bop 5 - 0? 
bottom :« i; 
repeat 

fiddle (top 4- bottom) div Bi 
if atmiddle] <« temp 
then top :« middle 
else bottom middle' 
uriil top + 1 « bo tt ft in 5 
for i s«* .f - 1 down to bottom do 
a C i + 11 s « a C i 1 5 
aCbottoml :« temp 

end 
end J 

C#***#4f***#^^***^** P Y-oc:edure shel lsort**********4M**iH*#^ 

procedure shel Isort (var a:list? n:inteqer>5 

var 

i> {index used for mass shiftinq of chain items* 

.is {some multiple of k* plus m index of item taei nq inserted} 

k» {becomes approx 1/3 'size of array * size of first increment) 
m, {.loop marker for each of the k chains} 
place* {.index if item in chain smaller than insertion spot} 
temp : inteqer? {item being inserted.! 
founds boolean; 
beqin 
k 1? 

while <3*k + 1 > < (n div 3) do 

k ■■ 3*k +lj {sets sise of increment based on n} 

repeat 

for a i" 1 to k do {for each of fche k chains} 
beqin 

H 8* k + mi 
while j <** n do 

beqin 

* temp i« aCj3; » element to be sorted} 

i J - k| 
found false? 



whi 3e not found and < i > 0) do 
if am > temp 

fchen 1 5- 3 -I: 

© 1 «se f nu nd 5 ss: t rue? 5 
p 1 ace : 138 i + k ; { where i tern is bei nq i ns&r tod * 
i : « i - k ? 
whiltt i >«■ place do 
beqin {move rest of chain down] 
aCi +k3 s* 8 aCil; 
i «w i - kj 
end ; 

atplace3 temp? 
i s ss i + k; 
©nd { whi le> 
end ; f for > 

k 8* (k - 1) div 3 {next increment si ;:e> 
until k f 1 
©nd i 

{*** ********** *****procedure quick«ort##**#*###******#***************> 
procedure qu icksort ( var as list; n: integer); 
procedure parti t ion <r*ss integer ;var j s inteqer >; 
var 

i s integer; 
temps inteqer; 
begin 

i' !* r + 1! finder one past pivot> 

i » 88 s; {last index in swap ranqeJ 

repeat 

while (aCU <« aCr3) and (i < s> do {find an element to tiwap} 
i l + 1 ; 

while <aCj3 >■ aCr3) and < ,j > r) do { *> 99 > 1 ,? } 

i i* i - U 
if i <J 

then 

beqin {swap} 

temp s 3 " aCil; 

a C i 3 j « a C j 3 ; 

aC.i] s** temp 
end; 

until i >« j5 {Swap finished for this pivot* 
temp 5= atr]; 
aCr3 s« aCj3; 

aC.il s~ temp; {puts pivot in place* {.} returned to callinq proqram 
end; 

procedure sort (m^ni inteqer ) ; 

{recursive procedure which redivides the sort field! 
var 

j 8 inteqer; 
temp s inteqer; 
beqin 

if n - m > 1 {provides the haltinq case* 

then 
beqin 

partition <m s n s j>i {receives the next dividing point, j> 
:RJC ,<5 * 



sort ai -i i - 1 > 8 
nor t < i i" t <n ) i 
end 

el&e if (n - m « 1 > and f*tm1 > aTnl) 
then 

begin *do a final swap) 

temp i« aCm3? 

aCml atn3? 

aCnl temp? 
end 

end.; 
begin 

sort < 1 y n) 
end 5 

* * * * * * # procedure-? bubble* *. •* * * * * * * * * * * •* * * * > 
procedure bnbh lesort < var as 1 ist 3ns integer ) ; 
var 

i> C index of items being compared > 

.is {loop variable limit} 

]*at«w) tch» fflaq for order in the data} 
temp s {swap temp} integer; 

beg i n 

j : « n - 1 ? 
repeat 

lastswitch :~ 1 ? 
f # or i 8«« 1 to .i do 
if atiJ > a Ci + 13 
then 

begin -Cawap} 

temp »« aCi3s 

ati3 :« ati * 13 J 

aCi 4- J 3 b« temp 5 

lastswitch :»« i $ 
end ; 

.1 i* lastswitch - 1 ? {end of unaorted data} 
until lastawitch » 1 

end ? 

{* * * # * * * * ^procedure selection* ********** 

procedure selection! var as 1 1st 5 n: integer ) ? 

var 

i 9 J ?n>k : integer S 
begin 

for j »« 1 to n - 1 do 
{generates atj3 thru aCn3} 
begin 
k *«* \ ; 

•fk carrirs smallest element in decreasing block) 
s 8« aC.i3; 

for i in j + l to n do 
begin 

if aH 3 < a then 
begin 

© i« aCi3j 
k 8« i i 
end J 



FRir 



52 



end ? 
«s t » a r 1 1 ; 
af.il 5** a ri: 1 5 
atkl 8 s38 © 
end ! 
end; 

{********* procedure heap* * ******** ******** } 
Procedure heap(var aslistSns integer) 5 
var 

in {index> 

ms integer ; {temp for array element;} 



procedure f i xheap ( var as list J top * hot torn 5 i nteqer ) 5 
{assume** i terns top + 1 to bottom ara a heap> 
{returns items top to bottom in a heap* 
var 

i ? { i nde* 3 

xainteqerji {temp for array & lenient > 

heqi n 

i a« & * top? 
if i <» bottom then 
beq in 

if i < bottom then if am <aCi + 11 then i i +1?{ largest child! 
if aCtop] < aCil then 
beqin 

x :« a C top 3 5 
artopn t» aCi3$ 

ati3 :« xj {heap is qood except at i> 
f ixheapfa* i 9 bottom) 
c*nd $ 

end; 
end ; 

beqinfheap} 

for i n div £ downto 1 do f i xheap<a» i ,n> J 
for i i« n downto S do 
beqin 

aC13s« aCi3; 
aCi3 :« xj 
fixheap(a«il!ii~l) 
end? 

end? 

{********* Main Program* ***** ************* y 
beq i n 

choice 05 

Writeln( 'This program sorts any number of items from & to ' v maxnu<n>| 

# Wr iteln( 'Enter the number to be sorted' >; 
readln(n) 5 

makelist<aiin) 5 

* showl ist (ajtn) i 
while choice <> 8 do 
begin 
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repeat 



t'Ji" i fc» 1 n ( 
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wr i t:P* 1 n ( 
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wr i teln( 


9 <+t 


b i nary i nser t ion* > ? 


writeln( 


p 5: 


Shell 'J $ 


writeln( 




selection 5 ' ) ; 


wr i teln( 


51 7: 


heap 51 ) ? 


wr i teln( 


> 8: 


end 11 ) ; 


readln( choice) 


• 

9 



until (choice >D) and (choice < 9) J 
case choice of 



t * beqin 

b 8 « * S 

bubb 1 e«snr 1 (b in) ? 
show! i * k ( h ? n ) ; 
end ? 
S: beqm 

b := a; 

quicksort (b * n) 5 
showl ist(b»n) ; 
end S 
3 n beqin 

b 8~ a? 

insert i on (b»r.) ; 
show list <b^n) 5 
end S 
'4 8 beqin 

b s« a; 

binaryinser t ion (b ? n) : 
showl ist <b«tn) ; 
end ; 
5t begin 

b s ~ a ; 

*v,ft* \ I sort (b $ n) 5 
showl ist (b *n) 
end 5 
6 8 begin 

b s« d! 

selection(b >n) ; 
showl ist<b >n> ? 
end i 
78 beqin 

b :« a; 
heap (bun) ; 
showl ist (bjiii) ? 
end ? 

end? 
* end 5 
end . 



